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Abstract 

A model for the dynamical evolution of a point charged particle interact- 
ing with higher order fields is described by an ordinary second order implicit 
differential equation. We show that such equation is free of run-away and pre- 
accelerated solutions of Dirac's type. The theory is Lorentz invariant, causal, 
and compatible with Newton's first law and energy-momentum conservation. 
However, it is not a local theory in the usual sense. In our model for point 
charged particles, covariant uniform acceleration motion is unstable and must 
be understood as an approximation. 

1 Introduction 

It is well-known the existence of fundamental problems in the electrodynamics of 
point charged particles interacting with its own radiation reaction field |15 [ 118 1 [21]. 
Fundamental investigations of the problem lead to the Lorentz-Dirac equation [8], 
that although based on basic assumptions, the Lorentz-Dirac equation is plagued of 
non-physical solutions. 

Some attempts to solve such problems range from reduction of order schemes 
for the Lorentz-Dirac equation [bij . re- normalization group analysis [201 121] . 
fluid dynamic models [3] , Feynman- Wheeler absorber theory of electrodynamics [S] , 
extended bodies theories [22j or higher order derivative field theories 0[T7], among 
other theories proposed. 

We think that the description of the radiation reaction of point charged particles 
must be considered from a new perspective. Our justification for seeking new ideas 
to solve the above problems relies on the fact that we think that the notion of field 
is physical meaningful (this excludes Feynman- Wheeler theory), that one should not 
base a fundamental theory like classical electrodynamics of point charged particles on 
a model like Lorentz-Dirac differential equation which is unacceptable because of the 
physical consequences that it implies (this excludes reduction of order theories) and 
on the idea that a point charged particle is a simple and nice model for elementary 
particles (despite the problems of infinite fields that this model has). 
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We seek for an alternative solution to the problem based on the simple model 
of a classical point charged particle and the fundamental principles of causality, 
covariance and energy-momentum conservation and locality is desirable. However, 
Dirac proved (and corroborated by many others) that under such hypothesis, the 
simplest dynamical equation must be the Lorentz-Dirac equation. Therefore, one 
needs to relax some of the assumptions or adopting a more complicated model of 
point charged particle, that still do not guarantee internal consistency. We choose 
to relax the notion of electromagnetic field and consider a theory with generalized 
higher order electromagnetic field^j, a notion that we also explain in this paper. 

The notion of generalized higher order field that we use can be justified as fol- 
lows. The standard definition of classical field as a section of a vector bundle over 
a manifold M presupposes its independence of the probe particle testing the field. 
However, we have learnt that such assumption can fail when dealing with atomic 
and subatomic systems, because those systems are fundamentally perturbed by the 
action of measurements. In a similar way, we can imagine that a notion of objective 
classical electromagnetic field (that is, a closed two form field living on the spacetime 
manifold M and independent of the probing particle) can fail to describe the effect 
of the dynamics of the probing charged particle. This happens if the effect of the ra- 
diation reaction is of a comparable magnitude than the electromagnetic field that we 
are measuring. Therefore, a notion of classical field that accommodates the existence 
of probing charged point particles is useful for a better understanding of feedback 
phenomena and in particular for radiation reaction in classical electrodynamics. 

In this paper we describe a differential equation for point charged particles that is 
free of the problems of the Lorentz-Dirac equation. Among other standard princi- 
ples, the derivation is also based in two new ideas: the notion of generalized higher 
order fields applied to the case of electromagnetic fields, and maximal acceleration 
geometry theory With the notion of generalized higher order fields on hand it 
is possible to eliminate the Schott term and other higher order than two derivatives 
from the equation of motion, using additional degrees of freedom from the higher 
order field. 

The notion of higher order field is not enough to provide an appropriated dynam- 
ical description of a point charged particle dynamics. However, we will show that 
using in addition using maximal acceleration geometries, one can obtain a model of 
generalized higher order electromagnetic fields compatible with the power radiation 
law of a point charged particle. Finally, imposing that the dynamical equation must 
be of second order, we obtain equation (|5.6p . Such equation is compatible with first 
Newton's law, it does not have un-physical solutions that plague the Lorentz-Dirac 
equation and is compatible with energy-momentum conservation. 

The procedure to obtain the equation (|5.6p breaks down in two situations. The 
first corresponds to trajectories of maximal acceleration. The second case corre- 
sponds to covariant uniform acceleration (that includes the first one but is qualita- 
tively different). In the case of covariant uniform acceleration, a similar analysis as 
in the non uniform case, led us to a differential equation which is a Lorentz force 
equation with a constant total electromagnetic field [10J. However, covariant uniform 
motion of point charged particles is unstable: in order to compensate the radiation 

2 For a more detailed mathematical description of the theory of generalized higher order fields, 
see the related work f3 . 
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reaction field and provide a total constant force, it is required an infinite precision 
in the calibration of the external fields. 

This paper explores the physical aspects of the equation (|5.6p as dynamical equa- 
tion for a point charged particle. The mathematical theory of generalized fields and 
differential forms, the fundamental elements of the geometry of maximal accelera- 
tion, as well as a mathematical theory of generalized higher order electrodynamics 
(not confuse with higher order derivatives theory like Bopp-Podolsky theory [2j [T7] 
are presented elsewhere. 

Notation. Through this paper M is a smooth manifold of dimension 4, TM is 
the tangent bundle, ./V := TM\{0} is the slit tangent bundle and r/ is the Minkowski 
metric with signature (—1, 1, 1, 1). Coordinates on M are denoted by (x° . X . X . X *j 
or simply by x^. We will identify points by its coordinates, since we are working on 
homeomorphic domains to R 4 . Natural coordinates on TM and N are denoted as 
(x,y). Further notation will be introduced along the paper. 



2 Notion of generalized higher order fields 
2.1 Jet bundles 

Jet theory is a natural framework for the study of the geometry of ordinary and 
partial differential equations (for a introduction to jet theory with applications in 
field theory see for instance [7J). Given a smooth curve x : I — > M, the set of 
derivatives (p = x(0), ^§)--->^§) determines a point in the space of the jet 
Jq(p) over the point p £ M, 

dx d k z 
da ' ' da' 

The jet bundle J fe (M) over M is the disjoint union 



Jo (P) : = { (z(0), ^lo, ^|o), VC fe - curve x : I -> M, x(0) = p € M, € / 



= u 


J k (x). 










dx | 


d 2 x | 


d k x | 






''do*' 



The projection map is 

\ : J k (M) -+ M, (*(0), ^| , ^| Q , ^| Q ) x(0). 

Given a curve x : I — > M, the fc-lift is the curve fc x : I — > Jq(M) such that the 
following diagram commutes, 

J k (M) 



I M. 



Each point on the lifted curve k x : I — > Jq(M) has local coordinates given by 
fc x : J-> J fe (M). 



(x(<t), ■••> ~~~s~tP~ s )-) where <r is the parameter of the curves x : I —> M and 
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2.2 Notion of generalized field 

In order to determine local properties (that is what one needs to write down a 
differential equation for point charged particles), we only need local expressions for 
the generalized higher order fields and to check that the equations are invariant 
under coordinate transformations. Let (M, 77) be the Minkowski space. 

Definition 2.1 The electromagnetic field F along the lift k x is a 2-form that in 
local natural coordinates is 

F( k x) = F(x,x,x,x\...) = (F f , u (x) + T fll/ (x,x,x,x,...))dx> 1 Adx v , (2.1) 
with F(x) G TA 2 M. The excitation tensor along the lift k x is the 2-form 

G( k x) = G(x,x,x, x, ...) = (G M „(x) + E^(x,x,x, x, ...))dx^ A dx v . (2.2) 
The density current in electrodynamics is represented by a 3-form 

J(x, x, x, x, ...) = (J^upix) + &^ U p(x, x, x, x, ...)) dx^ A dx u A dx p . (2.3) 

When acting on vector fields over M, generalized forms produce functions on J 7 ( Jq (M)). 
It can be shown that one can construct a Cartan calculus for such generalized forms 
in a similar way as the standard Cartan calculus of smooth differential forms. For 
instance, there is a coordinate free definition of the operator d^. The realization 
in local coordinates of the operator d^ is straightforward. For instance, if <fi is a 
generalized /c-form, its exterior derivative is 

d A 4> = d A {4>i(x(s), x, x, ...,x (fc) )) 
= d((f>i(x(s), x, x,...,x (k) )) 
= dj 4>i(x(s), x, x, x^) d^x-i A d^x 1 . 

With the aid of the nill-potent exterior derivative d^ and the Hodge star operator 
* defined by rj, one can write in complete analogy with the Maxwell's equations 
the following homogeneous and non-homogeneous equations. For the homogeneous 
equations one has 

d A F = 0. (2.4) 

Let us assume the simplest constitutive relation G = Then the inhomogeneous 
Maxwell's equations are 

d 4 * F = (p(J) + d 4 *T. (2.5) 

Using some constrains on T coming from compatibility with the equation of mo- 
tion of point charged particles, one can obtain an effective theory which is equivalent 
to Maxwell theory. The homogeneous equations are 

dF = 0. (2.6) 

Similarly, the non-homogeneous equations are 

d* F= J, (2.7) 
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and the conservation of the current density 

dJ = 0. (2.8) 

Therefore, one has an effective theory which is formally the standard Maxwell's 
theory. However, because the fields are originally in higher order spaces, there is 
still freedom to define a consistent particle dynamics. 

Generalized higher order fields can be physically interpreted as follows. The ideal 
way to measure an electromagnetic field is through its effects on probe point par- 
ticles. Usually the field is defined at the given point x G M, independently of the 
probing particle. This model idea works if the probe charged particle does not per- 
turbs considerably the field F{x). The situation changes if this is not the case, when 
the probe particle reaction field changes in sufficient amount the value of the field at 
the given point x. In this sense, the notion of a field living on M fails, since it is not 
physically observable. A generalized higher order field depends also on the probing 
particle trajectory, leading to fields F( k x) that depends on the probe particle curve. 
On the other hand, the flux of a field across ci surface element is measured by 
a standard macroscopic device. This is imprinted in our mathematical description 
by the fact that generalized higher order fields are forms acting on tangent vectors 
e /1 (x), e v (x)T TM defining the surface element A^ v . 



3 Maximal acceleration geometry 

Let D± be the covariant derivative in the direction x associated with the Levi-Civita 
connection of the metric r]. It induces a Levi-Civita connection 7^ vp on M and a 
non-linear connection on the bundle ir : TM \ NC — > M whose coefficients are given 
by 

N»„(x,y) = r vp y p , (3.1) 

with 7^ the Christoffel symbols and it : NC — > M the null-cone bundle. The 
inducing vertical forms are 

{5 y v = dy p + N 1 * v dx v , n,v = 0, 1, 2, 3}. 

The covariant acceleration vector field of the curve x : I —> M is given by D±x, 

(D ± x) m u := 2^ + 7^ up x u x p . (3.2) 

3.1 Elements of covariant maximal acceleration 

The covariant formalism for geometries of maximal acceleration that we use is de- 
veloped in [11]. The Sasaki-type metric on the bundle TM \ NC is the pseudo- 
Riemannian metric 

9S = r] pv dx p ® dx v + -3— r)^ Uy^ ® 5y u \ . (3.3) 

Let us consider a timelike curve x : I — >• M, parameterized such that for the lifted 
vectors g r 5( 1 i, 1 x) = —1. The metric (|3.3|) induces a bilinear, non-degenerate, sym- 
metric form g along the lift l x : I — > TM [llj . 
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Theorem 3.1 There is a bilinear, non-degenerate, symmetric form g along the lifted 
curve 1 x : I — > M, such for the corresponding horizontal lifted vector h x conformaly 
equivalent to rj, 

9^{x{t)) = f 1 — 2 J 77^. (3.4) 

■"max 

Proof. The tangent vector of the lift (x(r),x(r)) is (x,x). The metric gs acting on 
the vector field T = (x,x) at the point (x(t),x(t)) has the value 

g s (T, T) = (rj^ dx» ® dx u + -^—^ (V ® fy") ) [T, T) 

"■max 

1 

1 - -J2 ^ D ^ D * 3 

max 

which coincides with the value (g^ u dx 11 <S> dx u )(x, x), with g^ u given by (|3.4p . □ 
There is the following isometry, 

g s ( h X 1 , h X 2 )=g(X 1 ,X 2 ), (3.5) 

for any X 1 ,X 2 £ TTM. 

The generalized metric g will be called metric of maximal acceleration, because of 
the following property, 

Proposition 3.2 If g(x,x) < 0, then rj(D±x, D±x) < A^ax- 
Proof. From the expression ()3.4p and since g is non-degenerate, the factor 

, rj aX D x x°{t)D ± x x (t) . 

y 1 A2 ) 

"■max 

cannot be zero. Therefore, if g(x, x) < 0, then one has that r\{D± x, D± x) < A 2 max . 
□ 

Since the metric r/ is flat, there is a local coordinate system where the connection 
coefficients 7^ vp are zero (normal coordinate system) . In such coordinate system 
7^ vp = and therefore v = and the bilinear form (|3.4|) reduces to 

g{T) = (1 _ ^ d x» ® dx v . (3.6) 

■"-max 

In our considerations below we will use normal coordinate system associated with 
the metric r/. This is to avoid cluttering and because we are essentially disregarding 
gravitational effects. 



3.2 Perturbation scheme 

Let us consider a normal coordinate system associated with the Minkowski metric 
77. Let us denote by r the proper-time parameter along a given curve respect to g, 

t[x] = f y/-g(x'(t),x'(t))dt. (3.7) 
J t 
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We parameterize each world-line by the corresponding proper time r associated to 
g. The acceleration square function is defined to be 

a\r):= rj w x^x p . (3.8) 

Then the function e is defined in a normal coordinate system by the relation 

<r) = ( ti^ffl ) . (3.9) 

■"-max 

and in an arbitrary coordinate system by 

<r) : (3.10) 

mm 

The relation between 5 and rj is 

<?(r) = (l-e(r))r?. (3.11) 

e(r) determines a bookkeeping parameter eo by the relation e(r) = eo /i(r), where 

e = max{e(r), r € /}. (3.12) 

For compact curves, this definition always make sense. However, we will need to 
bound higher order derivatives in order to keep such parameter bound for non com- 
pact curves. Therefore, we can speak of asymptotic expansions on powers of O(e l ), 
with the basis for asymptotic expansions being {e l , I = —00, —1, 0, 1, +00}. 

Let us consider a normal coordinate system for n. We assume that the dynamics 
happens in a regime such that 

a\r) « A 2 max . (3.13) 

The curves X : I — >• M such that a 2 (r) = A%n ax are curves of maximal acceleration. 
We assumed that all the derivatives (e, e, e, ...) are also small. This is the way to 
stay away from the curves of maximal acceleration. 

The metric g defines different kinematical relations than n. If the parametrization 
of the world-line is such that g(x,x) = — 1 and the monomials in powers of the 
derivatives of e define a complete generator set for asymptotic expansions. 

In a maximal geometry space, the kinematical constrains are 

g(±,x) = -l, (3.14) 
g(x,x) = ^r](x,x), (3.15) 

g(xPx) + g(x,x) = + er,(x,x) (3.16) 

and similar conditions for higher derivatives obtained by derivation of the previous 
ones. Therefore, in a geometry of maximal acceleration (M,g), given the normal- 
ization g(x,x) = —1, the following approximate expressions hold: 



± p x p := g^x^xP = -l, (3.17) 
x p ± P := gp P x p x p =~+ O(e 2 ), (3.18) 
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4 A simple derivation of the Lorentz-Dirac force equa- 
tion 



4.1 Rohrlich deduction of Lorentz-Dirac equation 

Let us assume that the spacetime is the Minkowski spacetime (M, 77) . In a normal 
coordinate system of 77, the Lorentz-Dirac equation is the third order differential 
equation 

mx p = eF p v x v + \e 2 (x^ - (x p x p )x p ), := x p x a Vpa . (4.1) 

It has run-away and pre-accelerated solutions [8], both against what is observed in 
everyday experience and in contradiction with the first Newton's law of classical 
dynamics. 

We present a simple derivation of the Lorentz-Dirac equation (|4.ip . This deriva- 
tion is what we have called Rohrlich's argument [19] , It illustrates method Cartan's 
method [3], that we will use later in the contest of generalized higher order fields 
in combination with maximal acceleration geometry. One starts with the Lorentz 
force equation for a point particle interacting with an electromagnetic field F pu , 

m b x p = eF pv x\ (4.2) 

where mb is the bare mass and e the electric charge of the particle. Note that at this 
stage the bare mass is finite. Note that both sides are consistently orthogonal to x. In 
order to generalize the equation to take into account the radiation reaction, one can 
add to the right hand side of (|4.2|) a vector field Z along the curve x : R — > M. The 
orthogonality condition rj(Z(r), z(r)) = implies the following general expression 
for Z, 

Z p (t) = P p v (r){aix u {r) + a 2 x u (r) + a 3 x"), = 77^ + x p (t)x v (t), x p = 77^ 

(4.3) 

We can prescribe a\ = 0. Then using the kinematical relation x p x a r} pa = —x p x a r] p<J , 
one obtains the relation 

Z"(r) = a 2 x p (r) + a 3 (x p - (x p x^ pa ) x p )(r). 

The term a 2 x combines with the left hand side to renormalize the bare mass, 

(m b - a 2 )x p = mx p . (4.4) 

The argument from Rohrlich follows by realizing that in order to obtain the Lorentz- 
Dirac equation one needs 

a 3 = 2/3 e 2 . 

The same equation is obtained if instead of searching for a term containing the whole 
piece 



a 3 (x"(r) -(x p x° Vpcr )x»)(T), 



one requirers right hand side is compatible with the relativistic power radiation 
formula [El US], 

= - I eV*V)M *"(r). (4.5) 

In order to recover this relation, the minimal piece required in the equation of 
motion of a charged particle is —2/3e 2 (x p x p r] p(7 )x p . The Schott term §e 2 x is a 
total derivative. It does not contribute to the averaged power emission of radiation. 
However, in the above argument, the radiation reaction term and the Schott term 
are necessary, due to the kinematical constraints of the metric 77. Note that the 
above derivation is consistent with generalized fields with values on .F(Jq(M)). 

4.2 Rohrlich's derivation of the Lorentz-Dirac equation with max- 
imal acceleration 

Using the kinematical constraints for metrics of maximal acceleration, we can re- 
peat Rohrlich's argument. There is a small modification caused the bound in the 
acceleration. Repeating the procedure, one obtains 

Xl =l e 2 xPx p , A 2 = 0, A 3 = ^e 2 + O(e ), X k = 0, VA; > 4. (4.6) 

The corresponding modified Lorentz-Dirac equation is 

mx p = eF p v x u + \e 2 ("i"" - (x p x^)^) + O(e ), (4.7) 

with F^ „ := rj pp F p(J . This equation is formally identical to the Lorentz-Dirac 
equation. This fact proves that only maximal acceleration hypothesis is not enough 
to solve the problem of the Schott term in the Lorentz-Dirac equation. 

The same holds if we repeat the calculation in the framework of generalized higher 
order fields. Without the requirement of bounded acceleration one does not obtain 
a second order differential equation. In that case, one obtains again a differential 
equation that is formally the same than the Lorentz-Dirac equation. However, we 
will see in next section that if one combines maximal acceleration geometry with 
generalized higher order fields, it is possible to obtain a second order differential 
equation for point charged particles. 



5 A second order differential equation for point charged 
particles 

Let us assume that the physical trajectory of a point charged particle is a smooth 
curve of class C k such that g(x,x) = —1, x° > and such that the acceleration field 
is bounded from above. Using the generalized tensor fields one obtains the following 
general form for the differential equation, 

T^(x, x, x, x, ...) — B^Xi/ B u Xp -\- C^Xy Ci/Xp -\- F)pX v F) v Xp -)- ... x^ — Qnv& 
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Now we use the hypothesis that to first order on eo, Rohrlich's argument implies the 
Lorentz-Dirac equation^. This implies that we will have the expression 

m b x p = eF„+ (B^± u - x p B u ) x u 

+ (Cx u - x»C v )x u + {D*x v - x»D v )x v + 

and with F p v := g pp F pa = r] pp F pa . On the right hand side of the above expres- 
sion all the contractions that appear in expressions like {B p x v — x p B u ^x u , etc, are 
performed with the metric g instead of the Minkowski metric rj. The term 5m x p 
corresponds to the electrostatic mass due to the Coulomb field [15j.The other terms 
come from the higher order terms of the expression (|2. 1 j) of the electromagnetic field. 

The general form of a k-jet along a smooth curve x : R — > M implies the relations 

B"(s) = Px x"(s) + h x»{s) + (3 3 x p (s) + /3 4 ~x p (s) + • ■ •, 
C"( S ) = 7l i"(s) + 72 x»(s) + 73 x"( S ) + 74 + • • 

D"(a) = (Ji x"(a) + 5 2 x p {s) + 5 3 a: "(a) + 5 A x'^s) + ■■■. 

Let us assume that there are not derivatives higher than 2 in the differential equation 
of a point charged particle. One way to achieve this is to impose that all the 
coefficients for higher derivations equal to zero, /?& = 7& = 8}. = 0, k > 3. With 
this choice and using the kinetic relations for g, one obtains the expression 

m b x p = eF^ v x v - p 2 x fM + p 3 x" - -fak& - /3 3 (-a 2 (r)+ e')^. (5.1) 

We impose the condition that the differential equation governing the motion of a 
charge particle must be of second order and compatible with power radiation formula 
(|4.5|) . Therefore, if e(r) is different than zero, one obtains the relations 

fo= +^e 2 a 2 (s)j, (5.2) 
Pk = 0, VA; > 3. (5.3) 

Considering the leading terms in eo we obtain the differential equation for a charged 
particle in a higher order field, 

m b x» = eF p v x v \ e 2 a 2 (s) x" - \ e 2 a 2 (s) - x". (5.4) 

There is a re-normalization of the bare mass, 

n^b + - e a (s) - = m, e f 0. (5-5) 

Therefore, if e ^ we find the differential equation 

mi^eF,/- \ e 2 r tfXT x p x a x p , F p „ = g pv F pv . (5.6) 

The derivation of this differential equation is not valid when e = 0. We will consider 
this case separately. 

3 Strictly speaking, this not necessary. One can do the same analysis without mention Lorentz- 
Dirac equation and just follow Rohrlich's argument for the above general expression. 



10 



5.1 General properties of the equation (15. 6p 

Let us consider normal coordinates for rj and assume that we move on the domain 
of validity of the differential equation (|5.6p . Let us multiply equation (|5.6p by itself 
and contract using the metric g. Using the kinetic relations one obtains 

m 2 2 (1 _ £) = e 2 Flt p ±Pp u x -X (1 _ ^ + ( 2 e2)2 (a 2 )2 ^ ±v g ^ 

- 2eh 2 F\xPx u (l-e)r,^ 

= (1 - e) e 2 p W A i A ^ + (1 - e)(| e 2 ) 2 (a 2 ) 2 

- 2e^e 2 F%x p x u (l-e)7 ] ^ 

= (l- e )(^-(^e 2 ) 2 (a 2 ) 2 ). 
with the magnitude of the Lorentz force Fl given by 

F 2 L = e 2 F^F^x v x\ Xp . 



Proposition 5.1 For any curve solution of equation (|5.6|) one has the following 
consequences, 

1. The Lorentz force must be spacelike or zero, 

F 2 L >0. 

2. In the case the Lorentz force is zero, the magnitude of the acceleration is zero, 

Fl = ^ a 2 = 0. 

Proof. One can re-write the above expression as 

F 2 = ( 2 -e 2 ) 2 (a 2 ) 2 +m 2 a 2 , (5.7) 

from which follows the result. □ 
Therefore, one has the expression for u 



/ /In (^e 2 ) 2 

u=-m«C-i(l-\F L \J^ + c), C= [ -^J-. (5.8) 

Using this formula, one easily proves the following 

Proposition 5.2 The differential equation (|5.6p has the following properties: 

1. If there is not external electromagnetic field, the 4- acceleration is zero. 

2. If there is an external electromagnetic field, the acceleration is bounded by the 
strength of the corresponding Lorentz force. 
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3. For solutions of the equation (|5.6[) it holds the following asymptotic condition, 

lim F l (t) = lim a 2 = 0. (5.9) 

T — ^OO T— >0O 



Proof. From 1 to 3 are consequence of eq. (|5.8|) : the last consequence is by the 
hypothesis of maximal acceleration and the above properties. □ 

Note that a similar condition for the acceleration follows as soon as the external 
forces are zero F ext = 0. 

The differential equation (|5.6p is an implicit differential equation. One can solve 
x in the left side under some approximations. Equation (|5.6p can be read as 

2e 2 

x p (m5 p p + —Xpx^) = eF p v x v . 

In particular, if the radiation reaction term |e 2 a 2 x^ is small compared with the 
Lorentz force term one can treat it as a perturbation. As a first approximation one 
can consider x p = ~ F p x p and substitute in the equation ()5.6[) . The inversion of 
the operator 

M p p = (m5 p p + ^-x p F pX x x ) (5.10) 
is given by an expression in local coordinates as 

m 3m 2 

Therefore, the second derivative can be isolated 



°£= -^ K -^^F Kp x". (5.11) 



™ V F p v {^%- |J (5.12) 



This equation can be read as an approximation for (|5.6p . Once this has the ap- 
proximate isolated equation (|5,12p . one can use existence and uniqueness theorem 
of ODE theory [6] to state the following result, 

Proposition 5.3 Given two solutions z\ and Z2 of the equation of (|5.6p with the 
same initial conditions. Then they differ by a smooth function on the radiation 
reaction term |e 2 a 2 . 

Proof. Each of the solutions of equation (|5,6p can be approximated by a solution 
of equation of motion (|5.12p and the difference is given as power on |e 2 a 2 , starting 
at least at power 1. By uniqueness of solutions of ODE [6], given a fixed initial 
conditions (x p (tq) , x p (tq) is unique (in a finite interval of time r). Therefore, in 
such interval the two solutions differ by a polynomial on |e 2 a 2 . □ 

Note that equation (|5.12p can be written as a geodesic equation for a suitable 
connection coefficients. In a local holonomic frame {d p , g|jr}^ = o of TJq(M), 
let us consider a connection on TJq(M) with connection coefficients 

2 r^ va := ( 2 V^ d a y =^ va + K» va + V 1 va , (5.13) 
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with the tensors K and L defined by 

= -^(F^ VaX x> + F^ V(TX x>) * (5.14) 

1 1% 



L» V(T = -—F p(T FP v x». (5.15) 

All the rest of connection coefficients of 2 V are put equal to zero. Therefore, 

Proposition 5.4 The coefficients 2 V define a covariant derivative on Jg(M) and 
the geodesies 

2 V ± ± = (5.16) 
coincide with the solutions of the equation (|5.12|) . 

The geometric interpretation of equation (|5,12[) opens the possibility to investigate 
approximations to (|5.6[) using geometric and kinetic methods (see for instance [12]). 



5.2 Covariant uniform acceleration and mass renormalization 

The case of point charged particles such that e = requires special care. Our 
definition of covariant uniform acceleration is 

Definition 5.5 Covariant uniform acceleration curve is curve on the spacetime x : 
/ — > M such that along it the following constrain holds, 

d r](Dj.x, Dri) 

= 0. (5.17) 



Ht A 2 



This definition of covariant uniform motion is general covariant. Therefore, in order 
to make some computations we can adopt normal coordinates associated to rj. Then 
the kinetic relation (|3. 18|) holds. Therefore, although the four-acceleration is in 
general not orthogonal by g to the four velocity vector field, in the case of covariant 
uniform accelerated motion one has that g(x,x) = 0. By relation (|3.19p . 

g(x,x) = -g(x,x). 

Therefore, x is necessarily lightlike vector with the metric g, since the right hand 
side is necessarily zero, 

g(x,x) = (1 - e)rj(x,x) = (1 - e)e = 0. 

The implication of this fact is that the four-acceleration cannot be constant in a 
normal coordinate system, except for a 2 = 0. 

For maximal accelerations such that < oo our definition coincides with 

the standard definition of covariant uniform acceleration appearing in text books 
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[16] . It also implies the same notion of covariant uniform motion as solution of an 
homogeneous differential equation [ID] . 

Equation (|5.6p fails to describe covariant uniform motion. The reason for this is 
that in order to apply the mass renormalization procedure (|5.5p and the perturbation 
scheme, one requires that e ^ 0. Therefore, one needs to consider separately the 
case e = 0. Let us consider the following combination of equations (|5.5p and (|5.4p 



m b x^ = - eF* 1 v x" + fax* ■ 
Therefore, one has the following differential equation for covariant uniform motion, 

= ef v x v , m = m b {r) - /3 2 (r). (5.18) 

For covariant uniform acceleration one can see easily that 

d 

and this relation implies 



For trajectories with acceleration below the maximal acceleration A max , our defini- 
tion of covariant uniform motion of a point charged particle coincides with the one in 
[10] . In the case of r\ is the Minkowski metric and covariance is respect to Poincare's 
group. However, we will see below that this notion of uniform acceleration is highly 
unstable under perturbations. Note than at the times when tq e(r) = 0, tq must 
be an isolated critical point of a 2 . This means that one can extend the solutions of 
to such isolated points by continuity. 



5.3 Absence of run away and pre-accelerated solutions for the equa- 
tion 

Let us adopt a normal coordinate system associated to r\. We discuss the absence 
of non-physical solutions on (|5.6p . Run away solutions are solutions that have the 
following peculiar behavior: even if the external forces have a compact domain in 
the spacetime, the charged particle follows accelerating for ever. This is a patho- 
logical behavior of the Lorentz-Dirac equation. However, since the the asymptotic 
condition^] holds for equation (|5.6p . we shoe that our model is free of such problems. 

In order to investigate the existence of pre-accelerated solutions for the equation 
(|5.6|) , let us consider the example of a pulsed electric field [S] . For an external electric 
pulse E = (k5(t), 0, 0), the equation (|5.6p in the non-relativistic limit reduces to 

cix = k5{t), a=-^. (5.21) 
2e z 

The solution of this equation is the Heaviside 's function, 

ax =k, t > 0, (5.22) 
0, r < 0. (5.23) 
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This is not a pre-accelerated solution. In order to prove uniqueness on the space 
of smooth functions, let us consider two solutions to equation (|5,21|) . then they 
must differ by an affine function. This shows that in the non-relativistic limit the 
equation (|5.6p does not have pre-accelerated solutions of Dirac's type. Therefore, 
since the theory is Lorentz covariant, the equation ()5.6p does not have pre-accelerated 
solutions of Dirac's type. 

6 Discussion 

Combining a new notion of generalized higher order electromagnetic field with max- 
imal acceleration geometry, we have obtained the differential equation (|5.6p as a 
description of the dynamics of a point charged particle that takes into account the 
radiation reaction. Equation (|5.6p is free of run-away solutions and pre-accelerate 
solutions of Dirac's type. The hypothesis of maximal acceleration is necessary in 
order to keep under control the value of the acceleration. It also provides a book 
keeping parameter that allows to describe our perturbation theory. The hypothesis 
of generalized higher order fields is necessary too, since it provides the degree of 
freedom that we need in order to eliminate the Schott term in the Lorentz-Dirac 
equation. 

The use of generalized higher order fields for k = 3 and maximal acceleration 
above do not fix completely the dynamics of a point charged particle. However, the 
requirements that the dynamics is compatible with the covariant power radiation 
law (|4.5p . that the differential equation is second order, and the notion of field is 
extended minimally, completely fix the dynamics of point charged particles, except 
for covariant uniform motion. 

This different behavior of the uniform motion is because the method that we follow 
breaks down when e = 0. For covariant uniform acceleration one needs to consider a 
different differential equation ( in particular, following the same principles one arrives 
to (|5.18p ) than equation (|5.6|) . For a covariant uniform motion of a point charged 
particle, one needs to provide not only the initial conditions (x M (ro), x m (to)), but also 
to fix in anticipation the external field, such that along the world-line of the charged 
particle the total field is constant. Since the initial conditions (x^(tq), x^(tq)) can 
be used for both differential equations with in general different solutions, we should 
drop out one of them. 

The condition of covariant uniform motion requires of infinite precision on the 
determination of the external electromagnetic field in such a way that it compensates 
the radiation reaction field. Such a fine tuning is highly unstable, and the point 
charged particle will not behave uniformly accelerated under small fluctuations in 
the external field. Therefore, the stable description for the dynamics of a point 
charged particle is given by equation (|5.6p only. This argument and the assumption 
of requiring an stable rnode@ force us to abandon the notion of covariant uniform 
acceleration as a possible physical dynamics for point charged particles. Instead, 
it seems natural to consider the differential equation (|5.12p as an approximation to 
the motion of a point charge, when the radiation reaction is negligible compared 
with the external field. For instance, the difference on the solutions of (|5.12p and 

4 By a stable model we mean a weakly structurally stable pQ. 
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(|5.18p is given by the tensor IP" ua , which is by hypothesis in the approximation 
(|5.6|) = (|5.12p . a small perturbation. Therefore, for short time evolution and for 
initial conditions such that the radiation reaction is small compared with the total 
field Lorentz force, the difference on the solutions of the two equations is small and 
therefore, the difference between the solutions of ()5.6p and ()5.18p is also small. 

Note that the generalized fields are assumed to be well defined on the trajectory of 
the probing particle. This is contradictory with the fact that electromagnetic fields 
are infinite over a point charge particle trajectory, due to the Coulomb singularity. In 
order to consider this problem, we have adopted the method of mass renormalization. 
This method allows to consider well defined fields over any kind of trajectory, without 
the singularity of Coulomb type, which is renormalized. Therefore, all the fields and 
considerations in this model concern finite fields in the whole spacetime manifold 
M. 
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